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The nuclear many-body problem at the limits of stability is considered in the 
framework of the Continuum Shell Model that allows a unified description 
of intrinsic structure and reactions. Technical details behind the method are 
highlighted and practical applications combining the reaction and structure 
pictures are presented. 
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1. Introduction 

In this presentation we discuss specific features of the Continuum Shell 
O ' Model (CoSMo), the approach based on the projection formalism 1 , formu- 

lated in the classical book 2 and developed into a practical instrument in 
Refs. 3 ' 4 . The whole problem of many-body physics on the verge of stability 
has been extensively explored in the past, especially in relation to weakly 
bound nuclei. Alternative formulations and their first applications can be 



found, for example, in Refs. 5-7 . 

The goal of this paper is to highlight complimentary views on the nuclear 
many-body physics from the "inside" (structure) and "outside" (reactions) 
perspectives. The structure view is based on the traditional shell model 
where the effective Hamiltonian to be diagonalized plays the central role. 
New contributions to the effective Hamiltonian coming from the presence 
of continuum bring in non-Hcrmiticity and energy dependence. Overcom- 
ing these complications, it is possible to calculate in the same framework 
the cross sections of reactions, with their energy dependence and possible 
resonance behavior. The complementary picture that appears from the side 
of nuclear reactions is important for identifying resonances and compari- 
son with experiment. While the shell model approach to the many-body 
structure in discrete spectrum is firmly established the many-body reaction 
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physics is usually left for more phenomenological tools of the reaction prac- 
titioners. The purpose of Sees. 3.1 and 3.2 in this work is to accentuate on 
novel methods involved in calculation of Green's functions and associated 
time evolution operators that stay behind CoSMo. 

The example of a realistic application presented in the last section is a 
central point of the paper. The chain of helium isotopes is shown where a 
single picture combines different methods and different points of view on 
the same problem. The bound states of the conventional shell model below 
threshold are followed at higher energies by the solutions of the CoSMo 
effective Hamiltonian revealing resonances that coincide with the complex 
poles of the scattering matrix. The same resonances appear in the neu- 
tron scattering cross section plotted in the same figure. The discrepancies 
between the cross section peaks and resonance states emphasize subtle fea- 
tures of many-body dynamics in a marginally stable system. 

2. Structure 

Using the projection formalism one can eliminate the part of the Hilbert 
space related to particle(s) in continuum. This results in the effective Hamil- 
tonian TL that acts only in the "intrinsic" shell model space, 

H(E)=H + A(E)- l -W(E). (1) 

Here the full Hamiltonian H n is restricted to intrinsic space, and is sup- 
plemented with the Hermitian term A(E) that describes virtual particle 
excitations into excluded space and the imaginary term W(E) representing 
irreversible decays to the continuum. The new parts of the Hamiltonian 
(1) arc found in terms of the matrix elements of the full original Hamilto- 
nian that link the internal states |1) with the energy-labeled external states 
\c;E):A1(E) = (l\H \c;E), 

A 12 (£) = P.v. / dE'^ ^m^l , Wl2 (E) = 2n £ A^A?, 

c c(opcn) 

(2) 

Reduction of the effective space does not go without a price. The new 
properties of the effective Hamiltonian (1) are: 

1. For the description of unbound states the effective Hamiltonian is non- 
Hermitian which reflects the possible leak of probability from the internal 
system. 

2. The Hamiltonian has explicit energy dependence, making the internal 
dynamics highly non-linear. 
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3. The additional terms in the Hamiltonian that appear as a result of pro- 
jection can be complicated. Even with exclusively two-body forces in the 
full space, the many-body interactions appear in the projected effective 
Hamiltonian. 

By construction, the eigenvalue problem 



determines the internal part of the solution which is subject to the regular 
boundary condition inside matched to purely outgoing waves in the con- 
tinuum. For energies E below all thresholds, the amplitudes A\{E) vanish, 
and Eq. (3) determines discrete bound states with real £ — E. Above decay 
thresholds, Eq. (3) has no real energy solutions, and the stationary state 
boundary condition can not be satisfied. The similarity of this problem to 
that for the bound states makes it appealing to depart from the real axis 
and to find discrete non-Hermitian eigenvalues. The complex energy roots 
£ of (3) correspond to poles of the scattering matrix, see discussion below, 
and represent the many-body resonant Siegert states 8 . 

The transition into a complex energy plane may be rather impracti- 
cal, it causes computational complications related to numerous branch cuts 
and unphysical roots, the relation to observables becomes complicated and 
rather remote. As an alternative, the Breit-Wigner approach 9 is commonly 
used. Here the resonances £ = E — (i/2)T arc defined as Re [£ a (E)] = E 
and r Q = — 21m [£ a {E)\. In the limit of a small imaginary part (narrow res- 
onances) various definitions are equivalent. In the application of the CoSMo 
discussed below we use the Breit-Wigner approach. However, the general 
difficulty in parameterizing resonances in terms of centroid energies and 
widths should be noted. As demonstrated in Sec. 4, the problem becomes 
especially acute for broad resonances, high density of states, or in near- 
threshold situations. A look at the problem from the observable cross sec- 
tions is imperative. 

3. Reactions 

The picture where the nuclear system is probed from "outside" is given by 
the transition matrix defined within the general scattering theory 2 , 



The same transition amplitudes and propagation via intrinsic space drive 
the process shown schematically in Fig. 1. 



H{E)\a) =£\a) 



(3) 
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Fig. 1. Reaction process: the entrance channel 
b with amplitude A\ continues through inter- 
nal propagation started in the intrinsic state \2) 
driven by the non-Hermitian energy-dependent 
effective Hamiltonian (f) (with all excursions 
into continuum space included), and ends by 
exit from the intrinsic state |1) into the chan- 
nel a as described by the amplitude A^*(E). 

The poles of the transition matrix and related full scattering matrix 
S = 1 — 2iriT are the eigenvalues of Eq. (3) located in the lower part of the 
complex energy plane. The reaction theory is fully consistent with resonant 
description in Sec. 2. However, complexity of the many-body propagator in 
Eq. (4) with numerous poles, interfering paths and energy dependence can 
make the observable cross section which is a projection of poles onto a real 
energy axis quite different from a collection of individual resonance peaks. 

The solution of Eq. (3) with the large scale many-body Hamiltonian is 
a complicated task as extensively discussed in Refs. 4 . The calculation of 
the transition matrix (4) and of the cross section is yet another technical 
problem. The direct approach involving matrix inversion at all energies is 
extremely difficult and time consuming given large dimensions involved. 
The sharp resonances typically present in the spectrum make the process 
numerically unstable and require dense energy sampling to achieve a rea- 
sonable cross section curve. Absence of absolute numerical precision leads 
to instability near stable states embedded in the continuum where decays 
are prohibited by symmetry considerations. This problem is particularly 
troublesome within the m-scheme shell model approach. To overcome these 
difficulties, an alternative method has been developed which is discussed 
below. 



3.1. Unitarity and R-matrix 

The transition matrix (4) with the dimensionality equal to the number of 
open channels can be written as T = A^A, where the full effective Green's 
function Q (E) = 1/(E — H) includes the loss of probability into all decay 
channels. The factorized form of the non-Hermitian part W — 27tAAt in 
Eq. (1), where A represents a channel matrix (a set of columns of vectors 
A\ for each channel c) is the key for unitarity of the S matrix 10 . As 
shown in 11 , the simple iteration of the Dyson equation using the definitions 
Ti = H — iW/2 and G = (E — H)^ 1 leads to the following transition and 



a. 



b 



a 

Aj 



b 

A 2 / 



( E-H(E) ) 
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scattering matrices 



R S (5) 



1 + inR' 1 + inR 

where the matrix R = A^GA is analogous to the i?-matrix of standard 
reaction theory; it is based on the Hermitian part of the Hamiltonian H = 
H + A and computed as a function of energy using Chebyshev polynomial 
expansion in Sec. 3.2. 



3.2. Time evolution of the system and Green's function 

The technique behind the Green's function calculation in the CoSMo ex- 
tends the idea suggested in 12 where densities of states in molecular systems 
were computed using the Chebyshev polynomial expansion of the time- 
dependent evolution operator. 

First, the finite Hamiltonian matrix is rescaled and shifted by a constant 
so that the spectrum is mapped onto a generic energy interval [-1,1] using 
H — > (H — E + )/E-. The procedure involves scaling E_ and shifting E + 
parameters, E± = (E max ± _E m ; n )/2, that are determined by the upper and 
lower edges of the original spectrum E max and E m { n , respectively. Even 
within the traditional Lanczos diagonalization, the rescaling, although not 
required, is useful for providing numerical stability. Given a trivial nature of 
the rescaling procedure and its reversal, below we do not introduce special 
notations for the rescaled Hamiltonian. 

The energy representation of the retarded propagator is given by the 
usual Fourier image of the evolution operator, 

1 f°° 

G(E) = — - = -i dt exp(iEt)exp(-iHt), (6) 

E - H J 

where H is the Hamiltonian operator with a negative-definite infinitesimal 
imaginary part. The expansion factorizes the evolution operator using the 
Chebyshev polynomials as follows: 

oo 

exp(-iHt) = H) n (2 - S n0 ) J n (t) T n (H), (7) 

n=0 

where J n (t) is the usual Bessel function and the Chebyshev polynomials 
are defined as T„[cos(0)] = cos(n#). In comparison to the Taylor expansion 
or other methods evaluating the Green's function, the Chebyshev polyno- 
mials provide a complete set of orthogonal functions covering uniformly 
the interval [-1,1]. Although individual states can be resolved, the proce- 
dure is most effective when a significant energy region is involved, namely 
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for overlapping resonances. The asymptotic of the Bessel functions assures 
convergence of the series. The "angle addition" equations that follow from 
the definition of polynomials, 

(^)^m (x) — -^n+m(x) ~t~ ^n— m 

(x) , n > m, (8) 

are useful for the successive evaluation of series of vectors |A n ) = T n (H)\X) 
using the following iterative procedure: 

|A ) = |A), |Ai)=ZT|A), and |A n+ i) = 2H\X n ) - |A„_i). (9) 

In the CoSMo approach the calculations of reactions are performed using 
the Fast Fourier Transformation of Eq. (6), where the expectation value 
of the evolution operator in (7) is computed using iterative matrix-vector 
multiplications (9). In this way the R- matrix is computed which is then 
used to determine the cross section (5). The Chebyshev polynomials are 
divergent in the complex plane; therefore only the Hermitian part H = 
Hq + A corresponding to the R matrix can be used in the evolution operator 
(7). Using a conservative estimate it can be shown that n iterations lead 
to the energy resolution 4£L/n. Unlike the reorthogonalization problem 
in the Lanczos algorithm, lack of numerical precision in successive matrix 
vector multiplication does not leads to significant deterioration of the result. 
n = 1024 was typically used for CoSMo calculations. 



Fig. 2. Strength function of a 
dipole operator. Upper plot: effec- 
tive charges for protons and neu- 
trons are selected equal leading to 
a pure isospin T = CM opera- 
tor; only the CM states with en- 
ergies around 100 MeV have non- 
zero strength. Middle plot: effective 
charges are selected as -1 and +1 
for neutrons and protons, respec- 
tively, the resulting mixed opera- 
tor shows strength in both CM and 
non-CM states. Lower plot: the ef- 
fective charges are selected 
—N/A = -0.4 and e p = Z/A = 
0.6, which for 20 O excludes CM 
component from the dipole opera- 
tor; the resulting strength shows no 
CM excitation. 
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3.3. Center-of-mass separation 

To illustrate the effectiveness of the method, we make a digression from 
decays and continuum and discuss a stable large-scale shell model example 
of the center-of-mass (CM) problem. In Fig. 2 we show the strength function 
of the dipole operator. The strength function for a state |A) is defined as 

F X (E) = (X\6(E-H)\X) = --Im(A|G(£7)|A). (10) 

7T 

For the Hamiltonian H here, we consider the full s—p—sd—pf shell model 
space with positive parity states restricted to the sd shell, while the negative 
parity states include all one-particle-hole excitations from the sd shell. The 
two-body interaction is chosen as WBP 13 . A Lawson technique is used 
to address the CM problem with an artificial CM vibration Hamiltonian 
included into H with a large scaling factor. As a result, all states that 
correspond to the CM excitations appear at high energy, around 100 MeV 
in our example. The strength function of the dipole operator D = J2 a e a r a 
is considered, where e a is the effective charge of a particle a. In Fig. 2 the 
dipole strength for excitations from the + ground state of 20 O is plotted, 
namely Eq. (10) is evaluated with \Xd) — D\g.s). Depending on the choice 
of the effective charges for protons and neutrons, the operator D can be 
changed from the pure CM operator to the isovector operator containing 
no CM component. The change in strength of CM states is shown in Fig. 
2. 

4. Helium isotopes 

We conclude with a realistic example of CoSMo application to the chain 
of helium isotopes 4 He to 10 He that serves as an illustration of all tech- 
niques combined. The internal space of this simplified model contains two 
single-particle levels p 3 / 2 and pi/ 2 on top of the a-particle core. The sensi- 
tivity of the decay amplitudes A\ to the location of thresholds and to the 
parent-daughter structural relations lead to the necessity of considering the 
entire isotope chain. The effective shell model interactions were taken from 
14,15 xhese interactions are experimentally adjusted; thus it is assumed 
that the Hermitian renormalizations due to virtual particle excitations into 
continuum, A(E), Eq. (2), are already implicitly included; the energy de- 
pendence of A(E) is neglected. The diagonalization of the Hermitian many- 
body Hamiltonian within this valence space provides a conventional shell 
model solution. 
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One-body decays are accounted for in the model through the single- 
particle decay amplitudes defined as 

Al(E)= aj (e) (l;N\b]\a;N-l). (11) 

These amplitudes correspond to a single particle amplitude aj of the de- 
cay leaving a residual N — I nucleon state a, while the remaining nucleons 
can be seen as spectators. The amplitude function of energy is 

determined with the use of the Woods-Saxon potential that models the 
single-particle interaction between bound and continuum states. The pa- 
rameters of the potential are adjusted in order to adequately represent the 
4 He+n scattering. 

The two-body decays can be separated into sequential and direct ones. 
The sequential decays represent higher order processes generated by the 
same single-particle mechanism modeled here by the Woods-Saxon poten- 
tial. The direct decay requires introduction of new parameters describing 
instantaneous removal of an interacting pair. The model includes simplest 
two-body terms, see further discussion in 4 . 

In Fig. 3 the results of calculations are shown and compared to the 
experimental data 16-18 . The resonance states computed according to the 
Breit-Wigner definition are shown with discrete lines labeled with spin, 
parity and decay width. The same Fig. 3 contains a separate cross section 
calculation which implements techniques discussed in Sees. 3.1 and 3.2. The 
reaction calculation is performed with the use of the same Hamiltonian (1) 
and thus provides an important complementary picture to the resonant 
structure. The cross section shown is that for elastic neutron scattering off 
the ground state of the N — 1 nucleus. 

The model, in agreement with experiment, predicts ground states of 
4 ' 6,8 He to be particle-bound. The energies of these states by construction of 
the model exactly agree with the prediction of the traditional shell model. 
The states in the continuum arc approached from two perspectives, via 
the solution of Eq. (3) under the Breit-Wigner resonance condition and by 
directly plotting the cross section curve. The resonance centroids are shown 
with discrete lines with corresponding widths indicated. The continuum 
coupling changes the structure of internal states leading to resonant energies 
being in general different from those from the shell model prediction. The 
resonant patterns indirectly reveal information about structure of the states 
and dominant decay modes. The results for 7 He isotope agree with recent 
experiments 17 ' 18 . Our results support the "unusual structure" of the 5/2~ 
state identified by 18 . Due to its relatively high spin, this state, unlike the 
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Fig. 3. (Color online) CoSMo results for He isotopes. The states in the chain of isotopes 
starting from 4 He (top) up to 10 He (bottom) are shown as a function of the energy 
relative to 4 He. The horizontal dotted lines separate isotopes. For each case, the states 
from CoSMo are shown above experimentally observed states. The decay width (in units 
of MeV) along with spin and parity is shown for each state. The solid lines above CoSMo 
states show the elastic neutron scattering cross section from the spin polarized state of 
N — 1 isotope in the magnetic substate with M = (even) or M = 1/2 (odd mass) 
quantum number. 

neighboring l/2~ state, decays mainly to the 2 + excited state in 6 He. 

The comparison of the cross section curve with discrete resonances pro- 
vides a transparent picture revealing both usefulness and limitations of 
the resonant parameterization approach. The scattering cross sections start 
from thresholds set here by the ground state of the previous (N — 1) nu- 
cleus. The cross sections at sharp resonances, such as the ground state of 
7 He, agree well with the resonance parameterization. However, generally the 
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cross section curves are not symmetric and do not show a simple, Gaussian 
or Lorentzian, shape. The shape of low-lying states with widths big enough 
to reach threshold is particularly influenced. The association of the cross 
section peaks with the location of the resonance is ambiguous. A remark- 
able example is the case of l/2~ state in 7 He. The Breit-Wigner approach 
predicts an almost 3 MeV wide resonance at 2.3 MeV of excitation energy. 
The cross section curve, however, is only weakly influenced by such a deep 
pole and peaks near low energies reflecting primarily a proximity of thresh- 
old. This comparison of cross section and resonance parameterization may 
shed light onto the experimental controversy discussed in 17-20 . 
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